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Abstract 



We investigate an antiferromagnetic triangular Ising model with anisotropic 
ferromagnetic interactions between next-nearest neighbours, originally pro- 
posed by Kitatani and Oguchi (J. Phys. Soc. Japan 57, 1344 (1988)). The 
phase diagram as a function of temperature and the ratio between first- and 
second- neighbour interaction strengths is thoroughly examined. We search 
for a Kosterlitz-Thouless transition to a state with algebraic decay of correla- 
tions, calculating the correlation lengths on strips of width up to 15 sites by 
transfer-matrix methods. Phenomenological renormalization, conformal in- 
variance arguments, the Roomany-Wyld approximation and a direct analysis 
of the scaled mass gaps are used. Our results provide limited evidence that a 
Kosterlitz-Thouless phase is present. Alternative scenarios are discussed. 

PACS numbers: 75.30.Kz, 75.10.Hk, 05.50+q 
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I. INTRODUCTION 



Some spin systems display an intermediate critical phase, which is characterised by 
algebraic decay of correlations for a finite temperature range, and lies between the low- 
temperature ordered state and the high-temperature paramagnetic regime. This was first 
reported by Kosterlitz and Thouless (KT) [|]] in their study of the two-dimensional XY 
model, which exhibits long-range order only at T = 0. The KT transition from the paramag- 
netic to the critical phase is expected to involve an exponential divergence of the correlation 
length ||. 

Later it was found that some two-dimensional systems with discrete spin symmetry could 
accommodate critical phases as well. The first systems to be discovered in this category were 
the Z(N) models with iV > 4 Jose et al. @ found that, for the XY model with sixth- 
order anisotropy, the upper (Ti) and lower (T 2 ) limits of the KT phase are associated with the 
following values for the decay-of-correlations exponent T)(T): r/(Ti) = 1/4 and rj(T 2 ) = 1/9. 
This result has been used 0-§| as a criterion to establish the existence and limits of a KT 
phase in other systems, including the one which concerns us in the present work. 

Another class of spin systems with discrete symmetry displaying critical phases is ob- 
tained from models with a macroscopically degenerate ground state, by adding degeneracy- 
lifting fields or interactions. Such is the case of the nearest-neighbour Ising antiferromagnet 
on the triangular lattice. At T = this model has a finite entropy per spin 0, and al- 
gebraically decaying correlations ||. The behaviour of this system in a uniform magnetic 
field and with second-neighbour ferromagnetic interactions was studied by Nienhuis et al. 
||. They predicted that, when the nearest neighbour interaction, K nn , is set to — oo and the 
reduced field H = B/T is small, the algebraic phase persists, until it is destroyed at high 
fields. These predictions were verified subsequently (l^|Tl|. Other predictions, involving 



addition of finite (reduced) second-neighbour ferromagnetic interactions, K nnn , were made; 
it is expected that for large enough K nnn the algebraic phase exists also off the K nn = — oo 
line (i.e., for finite K nn ); when K nnn is increased further, however, a transition to a phase 
with long range order is expected. This is consistent with other results @|12[ (see below). 

The triangular Ising antiferromagnet with second-neighbour ferromagnetic bonds was 
initially proposed to model physical adsorption of gases on graphite |I3|| ; alternatively, it 
was used to describe the magnetic behaviour of CSC0CI3 and CsCoBr 3 JI4|.The model was 



predicted [12] to be in the universality class of the XY model with sixth-order anisotropy 
||, and thus to exhibit a KT phase. Monte Carlo results, together with the criterion of Jose 
et al., provided numerical evidence for this [|J. 

A phenomenological renormalisation group (PRG) analysis, however, produced clear 
signs of only one of the two transition temperatures associated with the KT phase JH) 



This discrepancy was ascribed to the small strip widths used, namely 3 and 6. In order 
to be able to analyse larger strip widths, Kitatani and Oguchi []5| introduced a slightly 
different model which incorporates only four of the six second-neighbour interactions per 
site. Each of the three second-neighbour sublattices then forms a square lattice, instead 
of the original triangular one. This introduces spatial anisotropy, with consequences which 
we shall discuss in what follows. In the original work of Ref. ||, the underlying hypothesis 
is that, since the two-dimensional character of the intrasublattice links was preserved, the 
qualitative behaviour ought to be in the same universality class of the model with all six 
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interactions per site. Those authors managed to study strips of width up to 12 sites with 
periodic boundary conditions (to which the Roomany-Wyld approximation [H| for the f3 
function |TJ] was applied), and 17 with helical boundary conditions (for which the spin- 
spin correlation functions were calculated directly and fitted to the criterion of Jose et al.). 
While results from the former approach were deemed inconclusive, the authors pointed out 
that the latter provided evidence of a KT phase. Subsequent Monte Carlo work was found 
to be consistent with this || by means of a finite-size scaling and renormalisation-group 
analysis of the simulation data, partly with help from the criterion of Jose et al.. Very 
recently, free energy and magnetisation were calculated for this model [IS] by a combination 
of transfer-metrix and mean-field ideas. The authors find an intermediate, incommensurate 
phase. The KT transition, in this interpretation, would be commensurate-incommensurate 
in character. An exact analysis of the ground-state of the model with second-neigbour bonds 
along only one direction also shows evidence of a critical phase ]19| ; when interactions along 



a second direction are included, so as to make the model identical with the present one, 
an exact treatment is no longer possible, even for the ground state. However, Monte-Carlo 
and transfer matrix data also point towards a commensurate-incommensurate transition at 
zero temperature {i.e. when the first- neighbour couplings are set at — oo) in this case, as 
the strength of second-neigbour bonds varies [J20| . The latter results are not to be directly 
compared to ours, as they pertain to the interacting domain-wall transitions which take 
place between ground-state configurations. 

Here we study the model proposed by Kitatani and Oguchi || on strips of width 3, 
6, 9, 12 and 15 with periodic boundary conditions. The interpretation of PRG results in 
searching for a KT phase is discussed, and these are shown to be of limited usefulness when 
considered on their own. The Roomany-Wyld (RW) approximation [16] for the f3 function 
[ I7J is implemented, for comparison with the corresponding results of Kitatani and Oguchi; 
a direct analysis of the ratio between the correlation lengths on strips of different widths 
(which, for quantum systems, correspond to the scaled mass gaps) [T^J2T|J2^] is given as well. 
We discuss the bearing of conformal invariance theory on the comparison of our results to 
previous ones. Finally we try to provide a consistent framework for our findings. 



II. PHASE DIAGRAM; SPECIFIC HEATS 

The model is characterised by the nearest-neighbour interaction J nn < and the second- 
neighbour coupling J nnn > 0, or alternatively by the ratios of these to the temperature, re- 
spectively K nn = J nn /T and K nnn = J nnn /T. A useful parameter is the ratio R = J nn / J nnn . 
Kitatani and Oguchi || restricted themselves to R = —2.0, whereas Miyashita et al || con- 
sidered both R = -2.0 and R = -5.0. We span the region -20.0 < R < 0.0. 

A phase diagram which, though not completely consistent with all our findings, appears 
to be quite plausible, is presented in Fig. [I]. Near the origin there is a high-temperature 
paramagnetic phase. Above the solid line we predict existence of a low-temperature phase 
with long-range order. The shaded region in between is our candidate for a critical (or 
massless) KT phase. Fixed- R lines, along which detailed studies were carried out, are also 
indicated. 

An examination of specific heats along lines of fixed R unveils a two-peak structure. This 
is especially apparent for R = —2.0, as shown in Fig. |2| for L = 3, ... ,15, for comparison 
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with Fig. 3 of Ref. |J, where L = 9 is the largest width used. However, when other 
values of R < —2.0 are investigated, the high-temperature peaks are displaced to ever 
higher temperatures and turn much broader, though those at the low end remain sharp and 
essentially unmoved. Fig. ^| illustrates this, with data for L = 15 and R = —2.0, —5.0 and 
—20.0 . The specific heat peaks apparently do not diverge with increasing L. Hence their 
positions do not identify transition points. The double-peak structure can be used only as 
a qualitative indicator of resemblance to Z(N) models with N > 4, which also exhibit such 
structure, and for which an intermediate massless phase is known to exist |2|]3|]. We shall 
not return to specific heats; instead, our work is focused mainly on correlation lengths and 
quantities derived therefrom, which will prove to be a safer ground for analysis. 

III. PHENOMENOLOGICAL RENORMALIZATION 

We begin by performing standard phenomenological renormalisation. The PRG pro- 
cedure consists of searching, for fixed R, for the fixed point {K* n , K* nn } of the implicit 
recursion relation 

£,L'(K nn , K nnn ) ^L(K nn , K nnn ) , 

L> ~ L ' {> 

where £l is the correlation length along a strip of width L. Usually best results are obtained 
for L'/L — > 1; as in the present case sublattice symmetry demands that these be multiples 
of 3, we use L' = L — 3. The correlation length is given by £l = 1/(C l°g A1/IA2I), where the 
geometric factor ( is 0,0] (2 /a/3) for K nn ^ and \/2 ioiK nn = 0; Ai^ are the two largest 
eigenvalues (in absolute value) of the 2 L x 2 L transfer matrix between adjacent columns of 
a strip of width L. As Ai is related to the strip free energy, it is always positive, while A2 
in the present case is one of a pair of complex conjugate eigenvalues with arguments ±27r/3 
(owing to the three-sublattice structure). If all six second-neighbour bonds were present 
||15|| , the transfer matrix would have to take into account adjacent pairs of columns, raising 
its dimension to 2 2L x 2 2L . Omission of two out of the six next-nearest neighbour couplings 
is thus technically worthwhile. One should, however, bear in mind the fact that, by doing 
this, the six-fold symmetry of the problem has been broken. Even though the model still has 
a six-fold degenerate ground state, the energies associated with domain walls between these 
phases depend now on the orientation of the walls, and one has to work out the symmetry 
of this new problem. It is quite clear that the criterion rj = 1/9 of Jose et al || to locate 
the low-temperature transition to a six-fold symmetric ordered phase should not be naively 
used here. 

For a system with a KT phase, one could expect Eq. |T] to exhibit a line of fixed 
points along a finite temperature interval, corresponding to the extent of the interme- 
diate phase P7| . Finite systems should display a precursor for such behaviour, such as 
having more than one solution to Eq. [I] at nearby points, in the region where a KT phase 
is expected. This does not happen in the present case, although the direct analysis of the 
scaled mass gaps (to be discussed below) shows that there are extensive temperature ranges 
along which other indicators of a massless phase occur. For each given R and L'/L, we 
find only one temperature as a fixed point of Eq. |l[ In Fig. |] we display the approximate 
critical lines given by the fixed points of Eq. [3] on the (K nn , K nnn ) plane, for the several 
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values of L'/L. Points on the extrapolated curve are obtained by keeping R constant and 
plotting the respective fixed-point couplings (e.g. K nnn ) for L'/L =6/9, 9/12 and 12/15 
against L~^, where ip is chosen to give the best straight-line fit through the three points, 
and is an estimate of the corrections-to-scaling exponent [p3| . We find that ip falls roughly 



in the range ~ 1-3 wherever such extrapolation is applicable (see below). 

We have also applied domain-wall renormalisation group ideas |24| by reversing the sign of 
the interactions that cross a line along the strip, which amounts to establishing antiperiodic 
boundary conditions f25| . The difference between the corresponding free energy and that 
for a strip with standard periodic boundary conditions gives the domain-wall free energy, 
which vanishes at criticality [F26J and is the quantity to be scaled. Again, only one fixed 
point was found for each R. The extrapolated phase diagram thus obtained (not shown 
here) is roughly the same as that coming from correlation-length scaling. Consequently, 
we did not pursue this line further. It is worth noting, however, that the finite-L data 
approach the extrapolation from the high-temperature side while those from standard PRG 
(see Fig. 01) come from below. This is to be expected from well-known duality relations 



between domain-wall energy and inverse correlation length [25] 



For isotropic systems, conformal invariance ||27] allows one to extract additional infor- 



mation from strip scaling, via the relationship between correlation-length amplitudes on a 
strip of width L at criticality and r/, the decay-of-correlations exponent : 

n = iM(r c ) . (2) 

Here, anisotropy is introduced by the missing second-neighbour bonds. Critical correla- 
tions are expected to decay with different prefactors in different directions, though with the 
same exponent. In some cases this can be explicitly dealt with p8|-|3"0|1, leading to a modified 



form for Eq. ^| which involves knowledge of correlation lengths along two different directions 
Z8fl . Unfortunately this it out of reach in the present model, whose very introduction was 



in order to enable one to build large- L transfer matrices along the specific direction where 
the bonds are missing. For K nn = or K nnn = 0, isotropy is restored and Eq. |] is valid, 
with results to be described below. 

Referring to Fig. EJ, the following comments are in order: 

(1) On the vertical axis (K nn = 0) one has three superimposed, disjoint, square lattices; 
thus the L arising in the conformal-invariance expression for r\ is one-third of the number of 
sites in the column. The corresponding critical quantities, extrapolated as above, are (exact 
values quoted in parenthesis): K nnn (0) = 0.441±0.001(0.4407...); y T = \jv = 1.00±0.01(1); 
rj = 0.250 ±0.001(1/4). 

(2) The extrapolated boundary leaves the vertical axis tangentially with a crossover ex- 
ponent ~ 1.55, where \K nn \ ~ (K nnn (0) — K nnn )^. This can be compared with the exact 
= 7/4 obtained by Slotte and Hemmer |31[ for the corresponding model with antiferro- 
magnetic second-neighbour interactions. 

(3) In the region between R ~ —0.01 and R ~ —0.3 the approximate critical lines become 
very close and cross each other twice. Thus, an extrapolation similar to that of Ref. [23 j is 
not feasible. 

(4) As R — > — oo the critical lines become horizontal; the extrapolated value is K nnn = 
0.313 ± 0.001. Already at R = —2.0 one is within less than 1% of that. Qualitatively, as 
\K nn \ grows the only ratio that matters is that between the two finite energies present, 
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namely J nnn and T. Such quantity is analogous to the magnetic field-to-temperature ratio 



in the zero-temperature nearest-neighbour triangular antiferromagnet [ 10| , |TT |. In order to 



check our results against those of Refs. [fT0| , jn]] , we set K nnn = and make \K nn \ large. 
This must coincide with their zero-field limit. Indeed we obtain t] = 0.4815, 0.4909, 0.4947, 
0.4965 respectively for L — 6, 9, 12 and 15, results in excellent visual agreement with Fig. 
2 of Ref. ]10[ and converging towards the exact value || rj — 1/2. 

The question remains of the interpretation of the critical line given by PRG in the context 
of KT transitions, where two phase changes are expected to occur. Following previous work 
|5|U, we try to fit our results to the criterion of Jose et al. 0. To this end we examine the 
exponent 77 on the critical line. Bearing in mind the above remarks concerning anisotropy, 
the quantity L/ir£,(T c ) does not in general correspond to 77; however, for large \K nn \ and 
finite K nnn one has essentially a triangular antiferromagnet with weak anisotropy induced 
by the square-lattice ferromagnetic bonds. In this region, corresponding to the horizontal 
section of the critical line in Fig. [|, Eq. |2| is expected to become asymptotically true. 
Fig. H shows L/7r£(T c ) as a function of K nn . It can be seen that for large |if nn |, L/ii^(T c ) 
becomes constant and extrapolates to 0.25 from below. On the basis of this, one would tend 
to state that the critical line of PRG marks the high temperature boundary T\ of the KT 
phase, at least in the region where the antiferromagnetic coupling is strong. This, however, 
is inconsistent with previous findings [|5],|6| (see below). Further, inspection of the 1 L vs. T 
curves for large \R\ shows no evidence of non-monotonic behaviour; in particular, no trend 
is found towards a second crossing close to temperatures where conformal invariance would 
give 7] near 1/9. An example is shown in Fig. |6|; for R = —5.0, at the couplings for which 
L/k^ = 1/9 we do not see any tendency for the different-L lines to cross, or even to approach 
one another. We then consider other quantities, in order to estimate bounds for a possible 
KT phase. 



IV. ROOMANY-WYLD APPROXIMATION 

To make contact with previous results, we investigate the Roomany-Wyld (RW) ap- 
proximation [|16j for the (3 function JIT], as done in Ref. || for R = —2.0. This is given 
by 

where £' denotes a derivative of £ with respect to temperature. Note that = at the 
fixed point of the PRG transformation (Eq. |l|), and one would expect it to remain at zero for 
an extended temperature interval corresponding to a KT phase. This would be in contrast 
with the monotonic behaviour exhibited in systems with an ordinary transition fi~7| . In 
Fig. we display f3f^ against temperature for R = —2.0 and L'/L=3/6, 6/9, 9/12 (already 
obtained in Ref. 0) and 12/15. Although the 3/6 curve is almost featureless, a trend arises 
towards flattening as L increases, especially when the 12/15 data are considered. The range 
of temperatures for which one would, from visual inspection, expect the L — > 00 curves to 
touch the axis, is in broad agreement with the approximate boundaries of the KT phase as 
obtained by previous authors 
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Fig. |8| shows Pll, against temperature for the largest strip width L = 15 [V = 12) 
and several values of R. It becomes apparent that the trend towards flattening noticed 
earlier turns even weaker, on either side of R = —2.0. As the RW approximation produces 
clear-cut plateaus in other cases where a critical phase is present (see e.g., Fig. 6 of ref. [0), 
one is led to concur with Kitatani and Oguchi's conclusion, by deeming the present results 
inconclusive for general R. 



V. SCALED MASS GAPS 

We then consider the ratio 

Q = (UT)/L)/(^(T)/L') . (4) 

This quantity was introduced in the study of the scaled mass gaps of quantum sys- 
tems [[Ll],[r[]], and was used in Ref. [|22| in the study of anisotropic Ising triangular anti- 
ferromagnets in an external field, which exhibit features similar to those under examination 
here. At the fixed point of PRG, Q = 1 and one expects it to remain close to this value for 
an extended temperature interval if a critical phase is present. This is numerically verified 
to a good extent in several cases |T6|j2l|j22| .In Figures |9| and |10| we show Q respectively for 



R = —2.0 and —5.0, the same values considered in Ref. M, and several pairs of strip widths. 



Figure [11] displays plots of Q against temperature for L — 15 and V = 12, and R = —0.5, 
-2.0, -5.0 and -20.0. 

Note that the fractional deviation of Q from 1 is much smaller than the corresponding 
shift of (3f^j from zero in Fig. As \R\ grows, the extent of the flattening region increases. 
For R = —20.0 and L — 15, Q rises quite steeply from 0.95 at T / J nnn = 2.4 to 1.0 at 
T / Jnnn — 2.96, the fixed point of PRG, and then remains below 1.005 up to T / J nnn = 10.0. 

The data from Q indicate that, at least for R < —2.0, there is a signature consistent with 
the presence of a KT phase at temperatures above that of the fixed point of PRG. This would 
support the results of Refs. |5|U . Visual extrapolation of the temperature range for which Q 
approaches 1 for large L would give the upper limits Ti/J nnn ~ 4.0 and ~ 5.2 respectively 
for R = —2.0 and —5.0. These compare well with the corresponding Monte Carlo results 
T\j Jnnn — 3.8 and ~ 5.0 . On the other hand, our evidence strongly suggests that the lower 
limit T2 must be the fixed point of PRG, whose extrapolated location is at T2/ J nnn — 3.19, 
both for R = —2.0 and —5.0 . This value is definitely higher than those quoted by Miyashita 
et al, resp. ~ 2.6 and ~ 2.75. 

For small \R\ the signature referred to above tends to disappear. If the KT phase is 
present at all, it must vanish as R — > 0, since at R = one has a simple structure of three 
independent, superimposed, ferromagnetic square lattices. In this context, one cannot rule 
out the possible existence of a multicritical point at some R t 7^ 0, such that for \R\ < \R t \ 
there would be only an ordinary transition. However, we have not investigated this aspect 
in detail. 

The trends arising from the examination of f3 RW and, especially, Q are difficult to recon- 
cile with the results for 77 coming from conformal invariance (though the latter are, admit- 
tedly, valid only for large \R\ when anisotropy becomes weak). One would expect 77 ~ 1/4 
at the high-temperature end of the massless phase, and i] < 1/4 at low temperatures still 
within the critical phase. This, however, is not seen in our results (recall, e.g. Fig. §). 
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VI. DISCUSSION AND CONCLUSIONS 



We have analysed four quantities for a wide range of values of the first-to-second- 
neighbour interaction ratio R in the Kitatani-Oguchi model: (a) the location of the fixed 
point of PRG; (b) the critical exponent rj, as given by conformal invariance; (c) the (3 func- 
tion in the Roomany-Wyld approximation; (d) the scaled mass gap Q. Overall, one can 
say that (c) does not provide conclusive evidence as regards the existence of a KT phase. 
On the other hand, the positions of the seemingly clear-cut plateaus arising from (d), above 
the fixed point of PRG given by (a) , do indicate existence of a massless phase; however the 
values of r\ obtained in this region, from (b), are inconsistent with the standard expectations. 

One possibility is that the model of Kitatani and Oguchi does not display a KT phase at 
all, owing to the fact (discussed above) that it is not six-fold symmetric. In this scenario, the 
transformation of the second-neighbour sublattices from triangular to square via the deletion 
of the appropriate bonds would be relevant in the renormalisation-group sense. For example, 
the work of Jose et al || indicates that, for the XY model with fourth-class anisotropy one 
has only one transition, with rj — 1/4. 

Although this is an appealing conjecture, we feel that our data from (d) for large \R\ 
definitely show a trend towards the onset of a critical phase. Based on this, we propose the 
approximate phase diagram shown in Fig. |T|, where the boundaries of the shaded area are 
obtained from extrapolation of the flat sections of the scaled mass gap diagrams. 

It is worth recalling that other generalisations of the triangular Ising antiferromagnet 
also display critical phases, without obeying the criterion of Jose et al. We refer specifically 
to the model with second-neighbour ferromagnetic bonds along only one lattice direction 
|T9| , and to one with continuous degrees of freedom superimposed on the original Ising ones 



32]. For the present model, we recall the zero-temperature transitions between domain 



states P0 |, and single out the Monte-Carlo renormalisation group results of Miyashita et al, 
summarised in Fig. 6 of Ref. ||, which do point towards the existence of a critical phase, 
and were obtained without recourse to the criterion of Jose et al. It would seem that this 
latter criterion is valid only for the model for which it was originally formulated and its 
closest relatives M with full six-fold symmetry. 
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FIGURES 



FIG. 1. Qualitative phase diagram. Heavy line is extrapolated from PRG (see Fig. 4). Shaded 
area gives massless phase, as obtained from scaled mass gaps. Constant-i? lines are shown, along 
which detailed studies were carried out. T marks (conjectured) multicritical point (see text). 

FIG. 2. Specific heat against temperature for R= —2.0. Bottom to top: L= 3, 6, 9, 12, 15 . 



FIG. 3. Specific heat against temperature for L=15 and R= —2.0, —5.0 and —20.0 . 

FIG. 4. Phase diagram in K nn -K nnn space, as given by the fixed points of PRG. The heavy 
line is an extrapolation from the 6/9, 9/12 and 12/15 curves. Some constant-i? lines are shown. 

FIG. 5. L/7t£l along the critical curves of PRG. The heavy line is an extrapolation from the 
6/9, 9/12 and 12/15 curves. 

FIG. 6. L/tt^l against temperature for R = —5.0 . Horizontal lines at 1/9 and 1/4 mark the 
values of rj at the boundaries of the KT phase for the XY model with sixth-order anisotropy. 

FIG. 7. Roomany-Wyld approximant fifjj against temperature for R = —2.0, and L/L'=3/6, 
6/9, 9/12 and 12/15 . 

FIG. 8. Roomany-Wyld approximant P^Jy against temperature for L/L'=12/15 and R= 
-0.5, -2.0, -5.0 and -20.0 . 

FIG. 9. Q = (£l (T) jV)j (£ L i (T) I V) against temperature for R= -2.0 and L/L'=3/6, 6/9, 
9/12 and 12/15. 

FIG. 10. Q = {i L {T)/L)/{i L ,{T)/V) against temperature for R= -5.0 and L/L'=3/6, 6/9, 
9/12 and 12/15 . 

FIG. 11. Q = (Z L (T)/L)/(fr(T)/L') against temperature for L/L'=12/15 and R= -0.5, 
-2.0, -5.0 and -20.0 . 
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